﻿ ACADEMiA ROMaNa iNSTiTUTUL NAtiONAL DE CERCETaRi ECONOMiCE iNSTiTUTUL DE PROGNOZa ECONOMiCa MODEL DE ESTiMARE A PiB‐ ULUi LUNAR Centrul de informare si Documentare Economica Bucuresti, 2005 Studiul de fata prezinta rezultatele partiale ale temei Analiza si prognoza evolutiei principalelor corelatii macroeconomice in Romania, comparativ cu alte tari din Europa din cadrul proiectului TiPOLOGiA CREsTERii ECONOMiCE; EVOLUtiA PRiNCiPALELOR CORELAtii MACROECONOMiCE iN PERiOADA DE TRANZitiE care face obiectului Contractului de finantare pentru proiecte prioritare PP4 S2 Programul National CERES incheiat de institutul de Prognoza Economica al Academiei Romane, cu Ministerul Educatiei si Cercetarii, prin institutul de Fizica Atomica in calitate de autoritate contractanta, la 25 Xi 2002 Durata programului noiembrie 2002-octombrie 2005 Subcontractant Academia de Studii Economice,Bucuresti Director de proiect : Marioara iORDAN Faza Pi iii 1 - Model de estimare a PiB-ului lunar Coordonator faza : Cristian Nicolae STaNiCa Cristian STaNiCa coordonator Lucian Liviu ALBU Mariana NiCOLAE Mihaela iONESCU Petre CARAiANi Carmen UZLaU Colaboratori externi Marin DUMiTRU Stelian STANCU Roxana CiUMARA Editat de CENTRUL DE iNFORMARE s i DOCUMENTARE ECONOMiCa REDACTOR-s EF - VALERiU iOAN FRANC SECRETAR GENERAL DE REDACt iE - AiDA SARCHiZiAN Redactor: ADELiNA BiGiCa Conceptie grafica, machetare si tehnoredactare: ViCTOR PREDA CiDE STUDii PROBLEME: Pro148 05 doc Redactia s i administrat ia: Bucuresti, Calea 13 Septembrie nr 13, sectorul 5, cod postal 76 117, telefon: 0040-1-411 60 75, telefax: 0040-1-411 54 86 Adresa postala: Bucures ti 5, casut a postala 5 - 72 Materialele cuprinse in acest buletin pot fi reproduse numai cu aprobarea conducerii institutului National de Cerceta ri Economice Volumele seriei pot fi identificate s i comandate fie in colectie anuala , respectiv iSSN 1222 - 5401, fie pe fiecare titlu in parte, respectiv pe iSBN alocat fiecarui volum Pentru volumul de fata: iSBN - 973 - 7940 - 85 - 7 CUPRiNS iNTRODUCERE 5 1 PREZENTAREA LiTERATURii PE PROBLEMA ESTiMaRii PiB-ULUi LUNAR 7 2 CONSTRUiREA UNUi MODEL TEORETiC DE ESTiMARE A PiB-ULUi LUNAR 9 3 ELABORAREA MODELULUi PENTRU CAZUL ROMaNiEi 18 BiBLiOGRAFiE 46 iNTRODUCERE Studiul "Model de estimare a PiB-ului lunar" investigheaza principalele metodologii utilizate in literatura pentru estimarea pe termen foarte scurt a indicatorilor macroeconomici de baza Doua dintre aceste metode (asa-numitele metodele algebrica si econometrica) sunt aplicate pe seria PiB-ului in cazul Romaniei din perioada 1997-2003 Estimarea PiB-ului lunar este utila in special pentru deciziile din sfera politicii monetare pe termen scurt in Marea Britanie, Comisia pentru Politici Monetare a recunoscut importanta PiB-ului lunar, inca de la lansarea acestui indicator economic in 1998, deoarece reprezinta unul dintre factorii determinanti ai ratei dobanzii Pe de alta parte, este recunoscuta relatia pe termen scurt intre modificarea ratei inflatiei si abaterea PiB-ului de la nivelul sau potential PiB-ul lunar vine sa completeze Conturile Nationale care sunt estimate cel mult trimestrial si cu o intarziere temporala La nivel lunar, datele disponibile (productia industriala) sunt insuficiente pentru caracterizarea evolutiei economiei "ca intreg" in consecinta, rezultatele analizelor economice pot fi eronate Avand informatii despre comportamentul output-ului economiei la frecvente inalte putem aplica modelul econometric la acest nivel pentru a prognoza indicatorii financiari si monetari De asemenea, informatiile continute in seria PiB-ului lunar pot fi utilizate pentru elaborarea unui model econometric care sa prognozeze PiB-ul lunar si trimestrial cu cel putin un trimestru inainte de aparitia datelor oficiale Lucrarea este structurata in trei capitole in primul capitol sunt prezentate pe scurt principalele teorii ale estimarii PiB-ului lunar din literatura de specialitate Originalitatea studiului consta in propunerea metodei algebrice pentru estimarea PiB-ului lunar in preturi curente, metoda ce vine in completarea metodelor deja cunoscute din literatura internationala Cel de-al doilea capitol detaliaza modelul teoretic al estimarii PiB-ului lunar si analizeaza trei metode: • Metoda algebrica de interpolare pentru serii trimestriale in preturi curente Avantajul metodei consta in evitarea calculelor econometrice restrictionate de testele statistice si utilizarea seriilor lunare si trimestriale in preturi curente mai usor accesibile • Metoda econometrica Chow s i Lin Avantajul metodei consta in simplitatea sa deoarece se bazeaza pe o ecuatie de regresie ai carei parametri sunt estimati cu ajutorul seriilor trimestriale in preturi constante Variabila dependenta este PiB-ul trimestrial Ecuatia cu parametri cunoscuti este aplicata ulterior la calculul PiB-ului lunar avand ca variabile explicative acelasi tip de serii accesibile la nivel lunar • Metoda filtrului Kalman de estimare a seriilor neobservabile Modelul teoretic al filtrului Kalman este prezentat pe larg in literatura internationala, dar aplicarea acestuia la cazuri particulare necesita stabilirea unor ipoteze de lucru specifice 6 Capitolul al treilea prezinta modelul de estimare si prognoza a PiB-ului lunar elaborat pentru cazul Romaniei Demersul stiintific este realizat in doua etape si anume: • Select ia s i prega tirea seriilor trimestriale s i lunare utilizate in model Analiza statistica a seriei trimestriale a PiB-ului in preturi constante, evidentierea problemelor metodologice de ordin statistic • Aplicarea metodelor algebrica si econometrica (Chow si Lin) pentru estimarea si prognoza PiB-ului lunar in economia romaneasca Compararea rezultatelor obtinute, elaborarea bazei de date 1 PREZENTAREA LiTERATURii PE PROBLEMA ESTiMaRii PiB-ULUi LUNAR Sensibilitatea sistemelor economice la socurile pe termen scurt din ultimii ani necesita tot mai mult elaborarea unor teorii care sa surprinda dinamica indicatorilor pe intervale scurte de timp Pentru aceste studii se impune imbunata- tirea bazei de date statistice Pe langa indicatorii trimestriali calculati in sistemul conturilor nationale exista cativa indicatori lunari insuficienti pentru a da o imagine completa asupra fluctuatiilor cererii si ofertei in acest sens, EUROSTAT in publicatiile sale recente si-a aratat interesul pentru estimarea lunara a principalilor indicatori macroeconomici (de exemplu produsul intern brut) pornind de la informatiile continute in seriile lor trimestriale Asa cum argumenteaza Lanning (1986), atunci cand utilizam serii cu date lipsa si dorim sa le completam sau dorim sa interpolam serii pentru a ajunge la o frecventa mai mare, avem la dispozitie doua metode de a rezolva aceasta problema O prima abordare consta in estimarea datelor lipsa simultan cu estimarea parametrilor modelului A doua posibilitate este construirea unei metode cu doua etape unde in prima etapa sunt generate datele lipsa, care devin astfel independente de algoritmul modelului in a doua etapa, seria obtinuta este folosita pentru estimarea econometrica a modelului Bazandu-se pe cercetari empirice, Lanning a aratat ca folosirea metodei in doua etape permite obtinerea unor parametri cu o varianta mai mica si deci prezinta o incredere mai mare fata de cazul estimarii modelului folosind seria incompleta Literatura legata de metoda interpolarii in doua etape poate fi clasificata in trei categorii (de fapt subiectul interpolarii este foarte vast iar realizarea unei clasificari obiective devine aproape imposibila) Prima contributie care se constituie ca o referinta importanta pentru toate abordarile ulterioare, apartine lui Chow si Lin si - care au prezentat o metodologie unica pentru interpolarea atat a seriilor de tip stoc cat si flux in prima etapa a metodei ei au estimat regresia seriei trimestriale (neobservabila la nivel lunar) in functie de variabilele trimestriale relevante care prezinta si serii de date lunare Apoi, in aceasta regresie cu parametri determinati au fost introduse seriile in valori lunare, ceea ce a permis estimarea seriei lunare neobservabile Aceasta metoda a fost folosita pentru interpolarea PiB-ului din Mexic de catre De Alba si a venitului populatiei din diferitele regiuni ale SUA de catre Schmidt Denton, Fernandez si Litterman au propus o abordare care minimizeaza o functie patratica a diferentelor dintre valorile seriei ce urmeaza a fi estimata si o combinatie liniara de serii observabile (numite serii aditionale) Aceasta metoda aproximeaza regresia lui Chow si Lin, dar se bazeaza pe ipoteze ceva mai complicate privind dinamica seriei care se interpoleaza precum si folosirea seriilor sub forma primei diferente O exemplificare pe date din Portugalia este data in Pinheiro si Coimbra 8 A treia abordare a fost realizata de Bernanke, Gertler s i Watson care au folosit un model al spatiului starilor bazat pe filtrul Kalman Avantajul metodei consta in faptul ca poate fi aplicata si la modelele neliniare Acest procedeu necesita stabilirea unor valori initiale pentru parametrii modelului precum si pentru seria neobservabila ce urmeaza a fi estimata O exemplificare pe date din Canada a fost realizata de Guay Originalitatea proiectului consta in propunerea metodei algebrice pentru estimarea PiB-ului lunar in preturi curente, metoda ce vine in completarea celor deja cunoscute din literatura internationala 2 CONSTRUiREA UNUi MODEL TEORETiC DE ESTiMARE A PiB-ULUi LUNAR Prima metoda, cea mai simpla dintre toate si mai aproximativa in ceea ce priveste acuratetea rezultatelor, a fost aplicata in cadrul activitatii de doctorat de cercetatorul Cristian Stanica ca o alternativa la metoda BNR de estimare a PiB-ului lunar (metoda BNR a stat la baza publicatiei Dihotomia real-nominal in economia romaneasca de tranzitie - Emilian Dobrescu ) Metoda directa se bazeaza pe ipoteza relatiei de proportionalitate dintre o componenta a PiB-ului si o serie reprezentativa din punct de vedere statistic, relatie valabila pentru lunile apartinand unui trimestru De exemplu, valoarea adaugata bruta din ramura industrie se determina la nivel lunar din productia industriala, coeficientul de proportionalitate fiind calculat din seriile trimestriale in acest studiu nu vom insista asupra metodei directe Metoda algebrica este un algoritm de calcul imbunatatit prin care determinam valorile seriei lunare necunoscute pornind de la o regula prestabilita de repartizare a valorilor lunare pe trimestru (seria este cunoscuta la nivel trimestrial) in cazul teoretic general, definim prin α indicatorul lunar neobservabil (VAB in industrie, VAB alte ramuri, sau PiB) si calculam valorile acestuia in interiorul trimestrului τ conform ecuatiilor: α3τ - 2 = α3(τ - 1) + d1,τ - prima luna α3τ - 1 = α3(τ - 1) + d1,τ + d2,τ - a doua luna α3τ = α3(τ - 1) + d1,τ + d2,τ + d3,τ - a treia luna (2 1) y3τ = α3τ + α3τ - 1 + α3τ - 2 - valoarea trimestriala a indicatorului τ = 1, 2,   T - T = numa rul de trimestre α3(τ - 1) se refera la ultima luna din trimestrul anterior iar d1,τ , d2,τ , d3,τ masoara variatia seriei neobservabile fata de luna precedenta Acesti parametri urmeaza a fi determinati din algoritmul de calcul in functie de ipotezele particulare impuse modelului in primul rand exista o relatie general valabila intre diferentele dintre trimestre notate cu ∆τ = y3τ - y3(τ - 1) = (notam) α(τ) - α(τ- 1) si diferentele lunare: ∆τ = 3d1,τ + 2d2,τ + d3,τ + d2,τ - 1 + 2d3,τ - 1 (2 2) Aceasta regula provine din conditia ca suma valorilor lunare sa fie egala cu valoarea trimestriala a indicatorului in al doilea rand, diferentele dk,τ se pot aproxima pe baza unor informatii suplimentare aduse de seriile aditionale la frecventa lunara, despre care teoria economica si statistica economica ne asigura ca sunt corelate cu seria lunara neobservabila Astfel, vom particulariza forma generala (2 2) a metodei algebrice considerand doua cazuri: cel al interpolarii cu serii aditionale, respectiv cel al interpolarii fara serii aditionale 10 interpolarea fara serii aditionale Rezultatele obtinute cu metoda algebrica fara serii aditionale se aplica la estimarea indicatorilor neobservabili (componentele PiB-ului) pentru care nu dispunem de nici o informatie ajutatoare la nivel lunar in ceea ce priveste rata de crestere ∆αt αt in acest caz, cea mai buna aproximatie se bazeaza pe ipoteza diferentelor constante, d1,τ = d2,τ = d3,τ = dτ , astfel incat am obtinut ∆τ = 6dτ + 3dτ - 1 sau, sub o alta forma, dτ = ( ∆τ   6) - (dτ - 1   2) (2 3) Este important sa facem precizarea ca diferentele dτ calculate cu ecuatia (2 3) pot fi si negative; cel mai adesea se intampla in trimestrul i pentru ca, de regula valorile PiB-ului in trimestrul iV sunt mai mari comparativ cu valorile din trimestrul i al anului urmator; in astfel de situatii valorile lunare ale indicatorului neobservabil vor scadea treptat una fata de cealalta in trimestrul i Relatia de recurenta (2 3) a pasului dτ permite determinarea acestuia numai daca se cunoaste valoarea lunara initiala d1 precum si diferentele trimestriale ∆i: d = 1 ∑ (-1)τ- i ∆i + (-1)τ- 1 d1 τ ≥ 2 (2 4) τ τ- iτ- 1 6 i = 2 22 Exponentul (τ - 1) de la numitor duce la convergenta rapida spre zero a parametrului d 1 2 2τ - 1 intr-o perioada de 1 an - astfel, indiferent de conditiile initiale, rezultatele modelului nu vor fi afectate De altfel, ne vom folosi de aceasta proprietate importanta a modelului luand pentru primul trimestru al seriei de date valorile lunare egale α1 = α 2 = α 3 = α(1) 3 (d1 = 0) Ecuatiile metodei algebrice in ipoteza diferentelor constante pentru trimestrul τ se simplifica considerabil; dupa cum urmeaza: τ3τ- 2 - 1 = - dτ τ(τ) 3 τ(τ) 3τ- 3τ- τ( τ ) 3 τ3τ- 1 1 = - = 0 τ( τ ) ponderile variaza simetric in jurul mediei (2 5) τ3τ 1 = dτ τ(τ ) 3 τ(τ ) iar dτ = ( ∆τ   6) - (dτ - 1   2) interpolarea cu serii aditionale Exista si situatii in care dispunem de serii lunare corelate cu seria care se interpoleaza Acesti indicatori lunari xt, care explica intr-un anumit procent (in functie de ponderea lor) evolutia indicatorului neobservabil α t, ii vom numi ca fiind "proxy" pentru α t - de exemplu, productia industriala PRODi este proxy pentru VABi (indicele ‘t' se refera la luna, t = 1, 2,  , 3T) 11 Pentru construirea algoritmului de interpolare cu serii aditionale vom considera ca αt este o functie de xt si o alta variabila independenta zt Atunci, la nivel lunar, unde variatiile variabilelor sunt mici comparativ cu valorile absolute, se poate aplica binecunoscuta regula Taylor de aproximare: ∆αt = ∂α∆xt + ∂α∆zt sau αt - αt - 1 = θt (xt - xt - 1) +  dt (2 6) t∂xt∂zt t t - 1 t t t - 1 t Astfel, metoda algebrica prezentata are avantajul ca tine cont si de influenta factorilor necunoscuti zt asupra modificarii variabilei neobservabile prin intermediul parametrului d t θt este panta functiei lui α t in raport cu xt cand ceilalti factori determinanti raman neschimbati Datorita lipsei informatiilor, cea mai buna aproximatie se obtine stabilind lips parametrii θt si d t constanti pe lunile unui trimestru, astfel regasim ecuatiile (2 1) in urmatoarea forma avantajoasa: α3 3 3 - 2 - 1 = θx⎛x 3τ-2- 1⎞ - d%τ α(τ) - 3 = θτα(τ)⎜⎝x( τ) - 3⎠⎟ - α(τ) α3τ-1 - 1 = θ x⎛x 3τ-1- 1⎞ τ 3τ-1 - = θ 3τ-1- α(τ) - 3 = θτα(τ)⎜⎝x( τ) - 3⎠⎟ ( τ) (2 7) α3τ - 1 = θ x⎛x 3τ- 1⎞ + d% τ α(τ ) - 3 = θτα(τ)⎜⎝x( τ) - 3⎠⎟ + α(τ) ∆τ = [θτ x(τ) - θτ - 1 x(τ- 1)] - 3x 3(τ - 1) (θτ - θτ - 1) + 6d%τ + 3 d% 3(τ - 1) dτ - 1 Totusi, acest model general nu poate fi aplicat decat in baza unor presupuneri specifice referitoare la valorile pe care le ia parametrul θτ in fiecare trimestru precum si la valoarea initiala d% (am vazut ca aceasta poate fi 1 considerata nula) O prima posibilitate ar fi sa alegem o valoare fixata θ0 valabila pentru toate trimestrele, pe care o determinam din regresia seriilor trimestriale α ( τ) s( τ) ( τ) i x Pe de alta parte, in situat iile in care indicatorul x nu este o componenta a indicatorului α(τ) in sistemul conturilor nationale, dar este variabila "cheie" care descrie evolutia lui α(τ), este indicat sa fixam panta θτ la valoarea α(τ) x(τ) caracteristica fiecarui trimestru τ Ecuatiile (2 7) se simplifica considerabil: α3τ- 2 = x 3τ- 2 - d% τ α(τ) = x( τ) - α(τ) α3 3 3τ- 1 = x 3τ- 1 α( τ) x( τ) α3 3 x 3τd% τ 3τ = x 3τ +dτ τ) ( τ) ( τ) (2 8) α(τ) x( τ) α(τ) d%τ 1x 3(τ- 1) ⎡⎛α( τ- 1) ⎞⎛x( τ) ⎞⎤d%τ-1 αd(ττ)= 21 x3((ττ-)1) ⎣⎢1 -⎝⎜α(τ) ⎟⎠⎝⎜xx(τ-1)⎟⎠⎦⎥ -2dατ(-τ1) 12 A treia posibilitate const a in alegerea θτ = 1 s i se aplica atunci cand (τ) (τ) indicatorul x este o componenta a lui α conform statisticii conturilor nationale (de ex sumele salariale brute sunt o componenta a valorii adaugate brute) Vom regasi in acest caz metoda diferentelor constante cu schimbarea de variabila αt   αt - xt :   [α3τ - 2 - x3τ - 2] - [α(τ) - x(τ)]  3 = - d τ [α3τ - 1 - x3τ - 1] - [α(τ) - x(τ)]  3 = 0 (2 9) [ α3τ - x3τ ] - [α(τ) - x(τ)]  3 = d τ dτ = [∆ (α(τ) - x(τ))   6] - (dτ - 1   2) Cele trei ipoteze ale modelului prezentate mai sus vor fi aplicate la interpolarea componentelor PiB-ului trimestrial in functie de natura acestora, iar prin agregarea componentelor se va obtine PiB-ul lunar Metoda Chow si Lin Chow si Lin au estimat PiB-ul lunar printr-o ecuatie de regresie avand aceiasi coeficienti ca si regresia PiB-ul trimestrial Variabilele explicative din regresia la nivel trimestrial sunt seriile agregate ale variabilelor lunare in forma matriciala, regresia seriei lunare neobservabile yt functie de cele m serii aditionale xt este y=Xβ+u, avand matricea de covariatie a erorilor V = E [uu′] iV iV iV La nivel trimestrial, se estimeaza parametrii regresiei y = X β + u Chow si Lin au gasit seria yˆ care aproximeaza cel mai bine PiB-ul lunar conform relatiei yˆ = Xβˆ + Λ(V) uˆ iV Seria are doua componente: cea obisnuita Xβˆ ( βˆ GLS GLS GLS este estimatorul GLS din regresia trimestriala) si o componenta Λ(V)uˆ iV de impartire pe luni a erorilor trimestriale uˆiV conform matricii de ponderare Λ (V) Aceasta din urma depinde de presupunerile facute asupra naturii matricii de covariatie V = E [u u ′] in modelul ales de noi consideram ca matricea de covariatie V a erorilor lunare este egala cu produsul dintre estimatorul dispersiei  ˆ 2 si matricea unitate i adica erorile au o distributie multinormala de medie zero Aceasta ipoteza conduce la simplificarea formei matricii Λ si la obtinerea erorilor in regresia cu serii lunare egale cu a treia parte din erorile estimate ale regresiei trimestriale Astfel, in acest caz particular avem toate elementele necesare pentru estimarea PiB-ului lunar pe baza informatiilor din regresia trimestriala Metoda filtrului Kalman consta in construirea unui sistem de doua ecuatii vectoriale care leaga variabilele lunare neobservabile α t (de ex PiB-ul lunar) de cele observabile zt (de ex o serie ipotetica lunara construita din date ale PiB-ului trimestrial) Prima ecuatie numita ecuatia de stare sau ecuatia de tranzitie descrie 13 (2 13) evolutia vectorului de stare α t continand cele n variabile neobservabile pe care vrem sa le estimam A doua ecuatie numita ecuatia de masurare leaga vectorul variabilelor observabile zt de vectorul de stare α t Considerand ca seriile contin T observatii lunare, t = 1, 2,  , T, modelul econometric are forma: α t = Ft αt - 1 + A t x t + R t u t (2 10) z t = H t α t + B t y t + N t v t (2 11) unde x t , y t , sunt vectori reprezentand variabile aditionale explicative (indicele productiei industriale lunare) care pot fi masurate, alese pentru a se imbunatati estimarea seriilor neobservabile Ecuatia de tranzitie este setata cu respectarea conditiilor initiale: α 1 = A 1 x 1 + R1 u 1 (2 12) iar vectorii erorilor sunt distribuiti multinormal: ⎛ut ⎞ ⎛⎛0⎞ ⎛Q 0 ⎞⎞ ⎜t⎟ " N ⎜⎜⎟ , ⎜ ⎟⎟ ⎝vt ⎠ ⎝⎝0⎠ ⎝0 G⎠⎠ unde Q si G reprezinta matricile de covariatie ale vectorilor erorilor in general, matricile parametrilor Ft, Ht, At, Bt, Q si G sunt estimate din conditia de maximizare a functiei de verosimilitate atas ata sistemului Calculele se fac iterativ pe masura ce se adauga noi date la seriile de timp, iar parametrii necunoscuti vor converge dupa un numar finit de pasi la valorile cautate Exista mai multe modalitati de estimare a modelului in functie de ipotezele impuse asupra erorilor Cele mai uzuale sunt ipotezele de normalitate (Harvey), dar se folosesc si alte procedee mai complicate cum ar fi: metoda proiectiei ortogonale (Brockwell, Chow, 1983), metoda generalizata a celor mai mici patrate (Sant) Aici vom aplica procedeele obisnuite bazate pe ipoteza normalitatii erorilor, iar in cazurile mai complicate cum ar fi regresia Chow, vom apela la metoda generalizata a celor mai mici patrate Pentru a estima PiB-ului lunar, in cazul economiei romanesti, vom construi un vector tridimensional αt = (yt, yt - 1, yt - 2 )′, in care variabila yt reprezinta PiB-ul lunar neobservabil Valorile PiB-ului trimestrial sunt introduse in variabila lunara z t dupa regula: z1 = 0, z2 = 0, z3 = prima valoare trimestriala a PiB-ului, z4 = 0 s a m d Matricile F, A, R se considera constante in acest caz, ecuatiile modelului se simplifica: α t = F αt - 1 + A x t + R u t (2 14) z t = H t αt + B t yt (2 15) cu conditia yt = x t + x t - 1 + x t - 2 Ecuatia de masurare este constrangerea ca suma valorilor PiB-ului pe lunile trimestrului sa fie egala cu valoarea trimestriala observata, ceea ce conduce la disparitia termenului eroare stocastic Pentru calculul seriei lunare a PiB-ului sunt supuse testarii doua tipuri de modele: cu serii aditionale (indicele productiei industriale) si fara serii aditionale 14 interpolarea fara serii aditionale Aceasta situatie este valabila numai daca exista suficienta informatie in seriile trimestriale astfel incat PiB-ul lunar sa urmeze un proces AR de ordin redus Tipul modelului AR este confirmat pe baza testelor cu serii trimestriale Procedura da rezultate pentru estimarea PiB-ului atat in preturi curente cat si in preturi comparabile Exemplu: Se presupune ca diferenta de ordinul i a PiB-ului lunar urmeaza un proces AR(1) stationar: ∆yt = φ∆yt - 1 + u t , u t fiind termenul eroare cu distributia N(0,   2) Daca exprimam seria in nivele, aceasta va descrie un proces AR(2): yt = (1 + φ)yt - 1 - φyt - 2 + u t , care se transcrie in reprezentarea spatiului starii: ⎜⎛yt⎟⎞⎛⎜1 + φ - φ 0⎟⎞⎜⎛yt-1⎟⎞ ⎜⎛1 0 0 ⎟⎞⎛⎜ut⎟⎞ ⎜yt-1⎟ = ⎜⎜ 1 0 0 ⎟⎟⎜yt-2⎟ + ⎜⎜0 0 0⎟⎟⎜ut-1⎟ (2 16) t-1 ⎜⎟t-2 ⎜⎟t-1 ⎝⎜⎜yt-2⎟⎠⎟⎜⎝ 0 1 0 ⎠⎟⎝⎜⎜yt-3⎠⎟⎟ ⎝⎜0 0 0⎠⎟⎜⎝⎜ut-2⎠⎟⎟ ⎝t-2⎠ ⎝t-3⎠ ⎝t-2⎠ cu αt = (yt , yt - 1 , yt - 2 )′ pentru t = 1, 2,  , T Ecuatia de masurare este conditia simpla de insumare a lunilor pe trimestre: z t = h t αt (2 17) unde h t = (0 0 0) pentru t = 1, 2, 4, 5, 7,  , T - 1, h t = (1 1 1) pentru t = 3, 6, 9,  , T Alt exemplu este sugerat de Bernanke, Watson (1997) Procedeul lor se aplica atunci cand PiB-ul (sau componente din PiB) este o serie i(1) cointegrata * cu o serie pt (poate fi un deflator) astfel incat seria lunara y t = yt  pt urmeaza un proces AR(1): y*t = φy*t - 1 + u t in reprezentarea spatiului starii modelul se transcrie: ⎜yt⎟⎛φ 0 0⎞⎜yt-1⎟ ⎛1 0 0⎞⎛⎜ut ⎞⎟ ⎜⎜y*t-1⎟⎟ = ⎜⎜ 1 0 0⎟⎟⎜⎜y*t-2⎟⎟ + ⎜⎜0 0 0⎟⎟⎜⎜ut-1⎟⎟ ⎝⎜⎜yy*tt 12⎟⎟⎠ ⎜⎝⎜ 10 01 00⎠⎟⎟⎜⎝⎜yy*tt 23⎟⎟⎠ ⎝⎜⎜00 00 00⎟⎠⎟⎜⎝⎜uutt 12⎟⎟⎠ z t = h t αt unde, comparativ cu exemplul anterior au fost redefinite: α t = (y*t , y*t - 1 , y*t - 2)′ pentru t = 1, 2,  , T si h t = (0 0 0) pentru t = 1, 2, 4, 5, 7,  , T - 1, h t = (pt pt - 1 pt - 2) pentru t = 3, 6, 9,  , T (2 18) (2 19) 15 interpolarea cu serii aditionale Metodele prezentate anterior pot fi imbunatatite prin adaugarea unor informatii suplimentare necesare interpolarii aduse de seriile aditionale Aceste modele sunt impartite in doua grupe: modele in care se considera ca PiB-ul lunar descrie un proces autoregresiv sau nu Exemplu: Chow si Lin au estimat PiB-ul lunar printr-o ecuatie de regresie avand aceeasi coeficienti ca si regresia PiB-ul trimestrial Variabilele explicative din regresia la nivel trimestrial sunt seriile agregate ale variabilelor lunare in forma matriciala, regresia seriei lunare neobservabile yt functie de cele m serii aditionale xt este y=Xβ+u, avand matricea de covariatie a erorilor V = E [uu′] iV iV iV La nivel trimestrial se estimeaza parametrii regresiei y = X β + u Chow si Lin au gasit seria yˆ care estimeaza cel mai bine PiB-ul lunar conform relatiei yˆ = XβˆGLS + Λ(V)uˆiV Seria are doua componente: cea obisnuita XβˆGLS GLS GLS ( βˆ este estimatorul GLS din regresia in serii trimestriale) si o componenta GLS Λ(V)uˆiV de redistribuire pe luni a erorilor trimestrialeuˆiV conform matricii de ponderare Λ(V) Aceasta din urma depinde de presupunerile facute asupra naturii matricii de covariatie V = E [u u ′] Modelul se transcrie in reprezentarea spatiului starii (Harvey si Pierse, 1984): α t z t = h t αt + b t (x t + x t - 1 + x t - 2 ) (2 21) in care s-a redefinit variabila de stare αt = y - x ' β si tt h t = (0 0 0) ; b t = 0 pentru t = 1, 2, 4, 5, 7,  , T - 1, h t = (1 1 1) ; b t = β′ pentru t = 3, 6, 9,  , T Rezultatele obtinute cu metoda regresiei Chow & Lin sunt aceleasi cu cele obtinute in formalismul spatiului starii Deoarece seriile xt au fost implicit incorporate in vectorul de stare αt, ecuatia de tranzitie este una de tip autoregresiv Metoda propriu zisa a filtrului Kalman este de maxima generalitate, avand avantajul de a face abstractie de orice presupunere asupra naturii matricii de covariatie V = E [u u ′] Forma ecuatiilor este simpla: ⎛ α ⎞⎛ yt- x'tβ⎞⎛0 0 0⎞⎛α ⎞⎛1 0 0⎞⎛ u ⎞ ⎜t⎟⎜tt ⎟⎜t - 1⎟ ⎜t⎟ = ⎜ αt-1 ⎟=⎜yt-1-x't-1 β⎟=⎜⎜1 0 0⎟⎟⎜αt-2 ⎟+⎜⎜000⎟⎟⎜ut-1⎟ (2 20) ⎝⎜⎜ αtt 21 ⎟⎠⎟⎜⎝⎜ytt 12-x'tt 12 β⎟⎠⎟⎜⎝⎜0 1 0⎠⎟⎟⎜⎝⎜αtt 23 ⎟⎠⎟⎜⎝⎜0 0 0⎠⎟⎟⎜⎜⎝utt 12⎠⎟⎟ 16 ⎛ x 1 ⎞ t ⎜⎛yt⎞⎟⎜⎛a1 a2 am ⎟⎞⎜ xt2 ⎟⎜⎛10 0⎟⎞⎛⎜ ut⎟⎞ ⎜yt-1⎟=⎜0 0 0 ⎟⎜t ⎟ +⎜⎜0 0 0⎟⎟⎜ut-1⎟ (2 22) ⎜⎜yt-1⎟⎟⎜⎜0 0 0 ⎟⎟⎜ M ⎟⎜⎜⎝0 0 0⎟⎟⎠⎜⎜ut-1⎟⎟ ⎜⎝yt-2⎟⎠⎝0 0 0 ⎠⎜⎜m ⎟⎟⎝0 0 0⎠⎝⎜ut-2⎟⎠ ⎜⎝ xt ⎟⎠ z t = h t αt (2 23) unde h t = (0 0 0) pentru t = 1, 2, 4, 5, 7,  , T - 1, h t = (1 1 1) pentru t = 3, 6, 9,  , T Termenul autoregresiv dispare din ecuatia de stare iar variabilele explicative apar in mod explicit Coeficientii ak se estimeaza in procedura filtrului Kalman din conditia de maximizare a functiei de verosimilitate Algoritmul Kalman de estimare a seriilor neobservabile Dupa ce au fost selectate modelele potrivite si au fost transcrise in reprezentarea spatiului starii, se trece la etapa estimarii parametrilor modelului si a seriei lunare neobservabile Considera m cazul general Filtrul Kalman este procedeul propriu-zis de calcul iterativ prin care se determina estimatorul eficient precum si matricea de covariatie a vectorului de stare α t conditionate de informatia disponibila la momentul (t-1) Daca αˆ t - 1 este estimatorul eficient (media conditionata) al lui αt - 1 pe baza observat iilor istorice pana la (z t - 1, x t - 1), atunci matricea de covariatie a erorii de estimare este conform definitiei din teoria statistica: Pt - 1 = E ⎡⎢⎣(αt-1 - αˆt- 1)(αt-1- αˆt- 1)'⎦⎥⎤ Pasii procedurii de iteratie Kalman sunt urmatorii: i La momentul (t-1) cunoscandu-se αˆt - 1 si Pt - 1 se fac predictii asupra estimatorului αˆ la momentul (t) pe baza ecuatiei de tranzitie: tt-1 αˆtt-1= Ft αˆt-1 + A t x t in timp ce matricea de covariatie a erorii de estimare se determina din relatia Ptt-1= FtPt - 1Ft' + RtQt Rt' tt-1 tt - 1ttt t Acestea sunt cunoscute sub numele de ecuatii de predictie ii Dupa ce au fost adaugate noile observatii z t la momentul (t) , estimatorul lui α t poate fi ajustat Ecuatiile de ajustare sunt: αˆ= αˆ + P Ht'Ot-1(zt - Htαˆ - Btyt) t⎪t - 1tt t tt⎪t - 1tt Pt = Pt⎪t - 1 - Pt⎪t - 1HtOtHtPt⎪t - 1 17 unde matricea O t este definita in continuare Estimatorul zˆ = Hαˆ + By poate fi interpretat ca media t⎪t - 1 ttt conditionata a variabilei z t la momentul (t) de catre observatiile istorice pana la (zt - 1, x t) Eroarea de predictie νt = zt - zˆ = Ht(αT - αˆ) + Ntvt ttt⎪t - 1 tT tt se numeste inovatie, deoarece reprezinta informatia noua cu care vine ultima observatie Vectorul ν t are media zero si matricea de covariatie Ot = HtP Ht' + Gt ttt⎪t - 1 tt iii Se fac in continuare predictii la momentul (t) pentru momentul (t + 1) si procesul se reia in forma finala a filtrului Kalman se reunesc ecuatiile anterioare intr-un singur sistem de recurente pentru a se face trecerea direct de la αˆ - la αˆ tt-1 t+1⎪t Pentru efectuarea calculelor se pot utiliza programele Eviews 4 0, Gauss, Stamp, Matlab Ecuatiile filtrului Kalman prezentate mai sus se bazeaza pe ipoteza ca perturbatiile si vectorul de stare la momentul initial (t = 1) sunt normal distribuite Aceasta permite calcularea de o maniera recursiva a distributiei vectorului de stare αt la momentul (t) conditionata de toata informatia existenta pana la momentul respectiv Distributia vectorului α t va fi ea insasi normala si astfel poate fi complet specificata de media si de matricea de covariatie calculate Daca nu mai este garantata ipoteza normalitatii, atunci nu se poate spune ca filtrul Kalman va calcula media conditionata pentru α t, dar va da in continuare un estimator eficient 3 ELABORAREA MODELULUi PENTRU CAZUL ROMaNiEi Selectia si pregatirea seriilor trimestriale si lunare utilizate in model Estimarea PiB-ului in preturi constante s-a facut pentru perioada 1997-2003 deoarece dispunem de o baza limitata de date trimestriale Toti indicatorii, atat cei de partea ofertei cat si cei de partea cererii, sunt exprimati in preturile medii ale anului 1997 Avand seriile in preturi curente, se obtin seriile in preturi constante prin deflatarea acestora cu indici de preturi estimati Mentionam ca datele 1 statistice pentru perioada 1997-2002 au fost calculate de drd Stanica Cristian in cadrul activitatii de doctorat sub coordonarea stiintifica a domnului academician Emilian Dobrescu si au dus la realizarea bazei de date trimestriale, versiunea iunie 2003, pentru Macromodelul Economiei Romanesti Pe aceasta cale, a fost realizata o colaborare stiintifica de cercetare deosebit de importanta cu specialistii Directiei de Conturi Nationale (Director Silvia Caragea; Director Gheorghe Mihai) din institutul National de Statistica Ulterior, in cadrul proiectului CERES actual, seriile au fost completate cu anul 2003 pe baza informatiilor furnizate de buletinul lunar al iNS Tabelul nr 1 Modificare procentuala din noul sistem fata de vechiul sistem - preturi curente 98Q1 98Q2 98Q3 98Q4 1998 PiB 1,7 0,3 0,7 0,4 0,7 Consum final menaje 1,7 0,2 0,6 -0,2 0,5 Consum final administratii 2,7 2,8 2,7 2,7 2,7 Formarea bruta de capital -0,3 -0,3 -0,3 -0,3 -0,3 - Formarea bruta de capital fix -0,3 -0,3 -0,3 -0,3 -0,3 - Modificarea stocurilor 0,0 0,0 0,0 0,1 0,0 Export net -1,7 2,7 0,5 -1,4 0,0 - Export -2,6 -2,7 -8,5 1,5 -2,9 - import -2,5 -1,3 -6,2 0,5 -2,2 Trimestrele anului 1998 in preturi curente sunt evaluate dupa ambele metodologii, SEC79 respectiv SEC95 neexistand diferente numerice semnificative (in tabelul prezentat mai sus exemplificam PiB-ul pe utilizari) Racordarea seriilor din perioadele 1997-1998 si 1999-2003 la preturi constante a fost posibila retinand in calcule anul 1998 din metodologia SEC79 1 Au fost consultate publicat iile: Dihotomia real-nominal in economia romaneasca de tranzit ie, Emilian Dobrescu (1997), CESTAT si buletinul lunar al iNS 19 Este cunoscut faptul ca institutul National de Statistica calculeaza indicii trimestriali ai PiB-ului fata de acelasi trimestru al anului anterior, incepand din anul 1998 Transformarea trimestrelor in preturile constante ale anului 1997 a fost realizata cu ajutorul indicilor de pret calculati pe baza informatiilor furnizate de indicii trimestrului curent fata de trimestrul corespunzator din anul anterior Suma trimestrelor in preturi 1997 se inchide pe cifra anuala publicata a indicelui de volum an an anterior Trebuie mentionat si faptul ca impunerea conditiei de insumare a valorilor trimestriale pe an in preturi constante duce la modificari ale indicilor de volum trimestru acelasi trimestru din anul anterior comparativ cu valorile publicate Stocurile sunt calculate rezidual pentru a se pastra egalitatea dintre valoarea PiB-ului pe productie si utilizari in preturi constante Astfel, erorile aparute din egalizarea PiB-ului in preturi constante pe productie si utilizari se reflecta in distorsionarea seriei in preturi constante a stocurilor si implicit a formarii brute de capital (iar formarea bruta de capital fix este serie nedistorsionata) Diferenta majora dintre metodologiile de evaluare ale PiB-ului SEC79 si SEC95 apare la structura consumului final in noul sistem SEC95 se disting doua concepte ale consumului final: 1 Cheltuiala pentru consum final - se refera la achizitionarea de bunuri si servicii de catre un anumit sector a) Cheltuiala pentru consum final a gospodariilor este formata din: Cumparari prin comertul cu amanuntul si servicii comerciale prestate populatiei plus transporturi, comunicatii; cumparari de pe piata taraneasca Chiriile imputate pe locuintele proprietarilor; productia pentru sine (consumul din productia proprie si consumul de bunuri procesate in gospodarii) Salariile in natura, avantaje in natura ale salariatilor, prestatiile sociale in natura (compensarea preturilor etc ) Vanzarile reziduale si corectie teritoriala b) Cheltuiala pentru consum final a administratiilor cuprinde: • Cheltuiala pentru consum individual a administrat iilor: bunuri s i servicii achizitionate de pe piata si furnizate fara prelucrari gospodariilor; bunuri si servicii nonpiata produse de administratia publica in domeniile: invatamant, sanatate, securitate sociala si actiuni sociale, sport, cultura • Cheltuiala pentru consum colectiv a administrat iei publice: servicii publice generale, aparare nationala, ordine publica, sanatate publica, protectia mediului, cercetare stiintifica, dezvoltarea infrastructurii si economiei 2 Consumul final efectiv - se refera la achizitionarea de bunuri si servicii de catre sectoare destinate consumului individual si colectiv Consumul final efectiv al administratiei publice corespunde numai consumului colectiv, consumul final efectiv al administratiilor private este nul, consumul final efectiv al gospodariilor este cheltuiala pentru consum a sectorului plus cheltuiala pentru consum individual a administratiilor Metodele de estimare a PiB-ului lunar sunt aplicate pe seriile "cheltuielilor de consum final" ale gospodariilor, administratiei publice si administratiei private in ceea ce priveste algoritmul de estimare a PiB-ul lunar in preturi constante, cu ajutorul seriilor lunare aditionale, trebuie sa tinem seama de doua 20 aspecte mai importante Primul se refera la modul de selectie al seriilor aditionale xt in primul rand, acestea trebuie sa fie corelate cu seria lunara neobservabila pe care o estimam (de ex o componenta a PiB-ului) si sa contina informatie relevanta pentru a se evita efectul erorilor induse care isi fac simtita prezenta in procedura de interpolare Un criteriu de selectie la indemana este urmatorul: seriile lunare potrivite pentru interpolare sunt cele care agregate la nivel trimestrial aproximeaza cel mai bine evolutia PiB-ului trimestrial Pe acesta ne-am si bazat atunci cand am calculat PiB-ul lunar Amemiya sugereaza o strategie bazata atat pe considerente economico-teoretice cat si pe evidenta statistica intuitia teoretica ne arata adesea care sunt seriile potrivite si care este forma functionala a modelului de care avem nevoie Pe langa aceasta trebuie sa avem o metoda statistica pentru a alege seria care da "cel mai bun" rezultat Cele doua criterii permit sa se faca alegerea finala a seriilor care vor fi folosite in modele Al doilea aspect se refera la numarul seriilor aditionale pe care le folosim in modelele econometrice Literatura de specialitate in domeniu arata ca seriile aditionale vin intotdeauna cu un cost al introducerii zgomotului (adica al erorilor din procesul de estimare) in seriile interpolate, daca sunt in numar mare Pe de alta parte, nu exista inca o metoda statistica care sa ne certifice in ce situatii este mai indicat sa aplicam modele cu serii aditionale (Chow si Lin) sau modele autoregresive (Bernanke, Gertler, Watson) de tipul Kalman Aplicarea metodelor pentru estimarea PiB-ului lunar in economia romaneasca Metodele directa si algebrica calculeaza (aproximeaza) PiB-ului lunar in preturi curente pornind de la semnificatia pe care o are acest indicator in Sistemul Conturilor Nationale Pentru estimarea PiB-ului pe latura ofertei dispunem de serii trimestriale in preturi curente din anul 1997 pana in anul 2003 in functie de datele disponibile descompunem PiB-ul astfel: PiB = VABi + VABA + VABO + NTAX, unde • VABi este valoarea ada ugata bruta in industrie; • VABA este valoarea ada ugata bruta in agricultura ; • VABO reprezinta valoarea ada ugata bruta din celelalte ramuri inclusiv PiSB; • NTAX sunt taxele indirecte nete Meritul metodelor directa si algebrica este ca incearca sa interpoleze trei din acesti patru indicatori (exceptie face VAB din ramura agriculturii pe care o calculam separat) si sa dea in final PiB-ul prin agregarea componentelor Pe latura utilizarilor descompunem PiB-ul in PiB = CONSM + CONSP + CONSG + FBC + EN, unde • CONSM, CONSP, CONSG reprezinta cheltuielile de consum ale menajelor, administratiei private si administratiei publice; • FBC este formarea bruta de capital; • EN este exportul net 21 Vom interpola numai consumurile si exportul net Astfel PiB-ul lunar este estimat pe latura ofertei iar prin diferenta se obtine seria lunara a formarii brute de capital Deoarece metoda directa este generalizata de metoda algebrica, prezentam in caseta nr 1 doar cateva particularitati ale acesteia: Caseta nr 1 Metoda directa 1 PiB pe latura ofertei 1 VAB lunara din ramura industrie variaza proportional cu productia industriala indicii de valoare lunari ai productiei industriale se obtin din indicii de volum ai productiei industriale si indicii de pret ai productiei industriale 2 Prin colaborarea cu specialisti de la institutul de Economie Agrara, au fost stabilite ponderile lunilor in trimestru pentru VAB din ramura agricultura 3 VAB lunara din celelalte ramuri & PiSB, variaza proportional cu fondul de salarii brute al ramurilor respective 4 impozitele nete lunare pe produs variaza proportional cu (TVA+Accize+Taxe vamale) 2 PiB pe latura utiliza rilor 1 Se calculeaza indicii de valoare ai vanzarilor de marfuri pe baza indicilor de volum si indicilor de pret ai cifrei de afaceri din comertul cu amanuntul De asemenea, se calculeaza indicii de valoare ai serviciilor de piata pe baza indicilor de volum si indicilor de pret ai cifrei de afaceri din sectorul serviciilor comerciale Se determina indicii de pret pentru comertul cu amanuntul, respectiv serviciile comerciale, pornind de la structura cosului de consum publicata in Buletinele lunare de preturi La vanzarile de marfuri se adauga produsul tutun iar din servicii se scot transportul, comunicatiile si utilitatile publice Cunoscandu-se valorile trimestriale raportate de iNS si indicii de valoare lunari se estimeaza lunar vanzarile de marfuri si serviciile comerciale Consumul final lunar al menajelor se calculeaza pe baza ponderii vanzarilor de marfuri si servicii de piata in consumul final stabilita la nivel trimestrial (ponderea variaza intre 50-70%) 2 Consumul final lunar al administratiei publice variaza proportional cu cheltuielile materiale & cheltuielile de personal din bugetul de stat+bugetele locale Se aproximeaza cheltuielile materiale & cheltuielile de personal cu cheltuielile totale (care nu contin cheltuielile privind datoria publica) din grupa: autoritati publice, aparare nationala, ordine publica, invatamant, sanatate, cercetare stiintifica Acestea detin in medie 88% din cheltuielile materiale si de personal totale, fiind neglijabile subventiile, transferurile si cheltuielile de capital 3 Consumul final lunar al administratiei private se considera constant in volum pentru lunile din cadrul unui trimestru, variatia sa fiind data de indicele preturilor de consum 4 Exportul net se determina ajustand valorile din balanta de plati la sumele pe trimestre Din PiB-ul lunar calculat prin metoda productiei, se obtine prin diferenta valoarea lunara a formarii brute de capital 22 Seriile lunare aditionale de partea ofertei utilizate de metoda algebrica pentru imbunatatirea rezultatelor estimarii PiB-ului lunar sunt: 1 PRODi - Productia industriala Au fost calculati luand ca baza de raportare media lunara a anului 1997 indicatorii "indicele productiei industriale - serie bruta" si "indicele preturilor productiei industriale" Notam aici faptul ca in perioada 1997-2003 s-a restrans aria de cuprindere a indicelui de preturi la piata interna indicatorul considerat de noi proxy pentru interpolarea VABi in preturi curente este "indicele de valoare al productiei industriale" in baza 1997 obtinut prin inmultirea celor doi indici 2 SALO - Alegem sumele salariale brute la nivelul economiei, exclusiv industria si agricultura, ca proxy pentru interpolarea VABO, fiind totodata un element de valoare adaugata Totusi, nu putem face o estimare lunara precisa pentru VABO deoarece procentul sumelor salariale calculat la nivel trimestrial variaza intre 20% - 35% si este variabil in ultimii patru ani, datorita influentei crescute a profitului asupra valorii adaugate brute 3 TVAA - (TVA + accize + taxe vamale) este considerata proxy pentru interpolarea impozitelor nete NTAX Figura nr 3 1 Evolutia productiei industriale in corelatie cu valoarea adaugata bruta din ramura industrie 250 200 150 100 50 0 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Productia industriala (linie punctata) si VAB din industrie (linie continua) in perioadele 1997 2000 si 2001 2003 Evolutiile la nivel trimestrial ale seriilor aditionale comparativ cu seriile pe care le interpoleaza sunt prezentate in figurile 3 1-3 3 Pentru a elimina efectele de scala induse de fenomenul inflatiei am optat pentru doua reprezentari grafice corespunzand perioadelor 1997-2000, respectiv 2001-2003 Toti indicatorii sunt raportati la media anului 1999 pentru prima perioada, respectiv la media anului 2001 pentru a doua perioada 23 Figura nr 3 2 Evolutia valorii adaugate brute din celelalte ramuri (fara agricultura) in corelatie cu sumele salariale brute 300 250 200 150 100 50 0 Figura nr 3 3 Evolutia indicatorului (TVA-accize-taxe vamale) in corelatie cu impozitele nete 250 200 150 100 50 0 Seriile lunare aditionale utilizate pentru estimarea PiB-ului lunar de partea utilizarilor sunt: 1 RS - Vanzarile de marfuri cu amanuntul, considerate proxy pentru consumul populatiei Datele statistice ale vanzarilor de marfuri sunt estimate si publicate de iNS numai pentru perioada 1994-1996 Se includ si vanzarile din comertul cu autovehicule, motociclete, carburanti Pentru perioada 1997-2003 nu exista date valorice disponibile Pentru a estima comportamentul acestui indicator la nivel lunar, am construit o serie tip indice de valoare cu baza media anului 1997, obtinuta prin inmultirea indicelui iPC cu indicii de volum ai cifrei de afaceri a intreprinderilor din comertul cu amanuntul Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 VAB din celelalte ramuri (linie continua) si Sumele salariale brute (linie punctata) in perioadele 1997 2000 si 2001 2003 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 24 2 Consumul final al administratiei publice este determinat in conturile nationale din cheltuielile salariale si materiale cuprinse in bugetul general consolidat la care se adauga cheltuielile de acelasi tip provenite din surse extrabugetare Cheltuielile bugetare aferente consumului sunt acoperite in proportie de 88% de cheltuielile totale (nu contin elemente de datorie, iar elementele de subventii si transferuri sunt mici) ale grupei autoritati publice, aparare nationala, ordine publica, invatamant, sanatate, cercetare stiintifica Le nota m cu GE s i le consider a m proxy pentru consumul guvernamental 3 Exporturile si importurile de bunuri si servicii din balanta de plati se ajusteaza direct la sumele pe trimestre, nemaifiind necesara aplicarea unei metode de interpolare De asemenea, si datorita faptului ca este neglijabil comparativ cu PiB, consumul administratiei private CONSP este considerat constant pe lunile unui trimestru, fiind egal cu media lunara Corelatia dintre variabilele CONSG GE si CONSM RS este indicata in figurile 3,4-3,5 Seriile CONSM RS au fost raportate la media anului 1999 pentru intervalul 1997-2000, respectiv la media anului 2001 pentru intervalul 2001-2003 Figura nr 3 4 Evolutia consumului guvernamental in corelatie cu cheltuielile guvernamentale ale grupei reprezentative 300 250 200 150 100 50 0 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Figura nr 3 5 Evolutia consumului populatiei in corelatie cu cifra de afaceri din comertul cu amanuntul 250 200 150 100 50 0 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Q1 Q2 Q3 Q4 Chelt Guvernamentale totale ale grupei reprezentative (linie punctata) si Cons Guvernamental (linie continua) in perioadele 1997 2000 si 2001 2003 25 Din analiza graficelor prezentate se desprind urmatoarele concluzii: • Productia industriala prezinta o evolutie corelata cu VAB din ramura industrie, cu toate ca in ultimii ani VABi are o variabilitate mai accentuata O corelatie semnificativa exista si intre consumul populatiei si cifra de afaceri din comertul cu amanuntul, cu abateri in primul si ultimul trimestru al fiecarui an Din punct de vedere al teoriei economice si al metodelor statistice aceste corelatii sunt justificate • Aproxima ri mai put in exacte se obt in in celelalte cazuri: cheltuielile guvernamentale ale grupei reprezentative urmaresc evolutia consumului guvernamental dar sunt mai volatile, iar in ultimii ani exista o discordanta semnificativa datorita caracterului provizoriu al conturilor nationale Pe de alta parte, sumele salariale brute nu explica abaterile mari din primul si ultimul trimestru al fiecarui an inregistrate de VAB din celelalte ramuri (exclusiv agricultura, industrie), datorita influentei crescute a excedentului brut de exploatare • Suma TVA+accize+taxe vamale aproximeaza satisfacator seria taxelor nete cu exceptia anului 2003; aceasta discordanta se datoreaza caracterului provizoriu al conturilor nationale care se observa si la nivel de an Caseta nr 2 Din observatiile anterioare ar rezulta faptul ca interpolarea PiB-ului pe luni are un caracter mai putin exact Totusi, conditia ca suma PiB-ului pe lunile unui trimestru sa fie egala cu valoarea trimestriala cunoscuta corecteaza erorile de interpolare iar introducerea seriilor aditionale in metoda de interpolare da un caracter specific evolutiei lunare a PiB-ului Consideram ca la nivel lunar factorii accidentali care abat seria PiB-ului de la modelul evolutiv specificat de seriile aditionale nu induc variatii bruste Factorii lunari accidentali sunt factorii necunoscuti; amintim aici factorii financiari care determina evolutia profitului pe luni; factorii agricoli care influenteaza marimea autoconsumului la nivel lunar; achizitiile guvernamentale si cheltuielile salariale din surse extrabugetare ale sectoarelor aparare si ordine publica Se poate mentiona ca estimarea PiB-ului la nivel lunar se face in jurul unui trend lunar potential determinat de maximul de informatie disponibila continuta in seriile aditionale Estimarea PiB-ului lunar in preturi curente se realizeaza cu ajutorul a trei tipuri de metode deduse din forma generala a metodei algebrice Astfel: Metoda algebrica directa cu variabile aditionale considera ca VABi, VABO, NTAX, CONSM, CONSG variaza direct proportional cu indicatorii proxy corespunzatori conform ecuatiilor (2 8), in care diferentele d τ sunt luate zero in fiecare luna Metoda pune accentul pe influenta factorilor cunoscuti ca si cum acestia ar determina singuri evolutia variabilei neobservabile La polul opus se afla metoda algebrica directa fara variabile aditionale in care modificarile lunare ale seriei interpolate din cadrul unui trimestru sunt constante, egale cu parametrul d%τ Estimarea lui d%τ se face pe baza ecuatiilor (2 5) luand ca unica informatie ajutatoare cunoasterea diferentelor dintre trimestre 26 indicatorul care induce cea mai mare incertitudine la estimarea PiB-ului lunar este VABA - valoarea adaugata bruta din agricultura Seria lunara a fost dedusa direct din valorile trimestriale cu ajutorul unor ponderi exogene stabilite pentru cele 12 luni ale anului, in urma discutiilor avute cu specialistii de la institutul de Economie Agrara Metoda algebrica generala tine cont atat de influenta factorilor cunoscuti, cat si de influenta celorlalti factori, conform ecuatiilor (2 8) Cuantificarea factorilor necunoscuti pentru lunile aceluiasi trimestru se face prin intermediul diferentelor constante d%τ Nu exista inca un raspuns definitiv la problema stabilirii valorii initiale pentru d% τ (valoarea de start a seriei) Un indicator care evolueaza in mod continuu, nu in salturi cum se intampla in economia romaneasca (apar diferente mari intre ultimul trimestru al unui an si primul trimestru al anului urmator) este foarte bine estimat daca a fost luat d% τ = 0 in primul trimestru al seriei de date Pare mai realist in cazul economiei romanesti sa interpolam elementele din PiB luand d% τ = 0 in primul trimestru al fiecarui an in lucrare sunt analizate ambele situatii Metoda generala este mai indicata deoarece gaseste un compromis intre influenta factorilor cunoscuti descrisa de seriile aditionale la nivel lunar si influenta celorlalti factori necunoscuti descrisa de variatiile indicatorului la nivel trimestrial Comparativ cu metoda algebrica directa, seriile VABi, NTAX, CONSM, CONSG sunt "imbunatatite" conform modelului (2 8), iar pentru seria VABO se aplica modelul (2 9), deoarece SALO este element de valoare adaugata Seriile VABA, CONSP, EN raman nemodificate PiB-ul lunar se determina prin insumarea elementelor de pe latura ofertei iar prin diferenta se obtine seria lunara a forma rii brute de capital FBC PiB-ul lunar estimat prin cele trei modele este prezentat in figura nr 3 6 Pentru a putea fi vizualizate mai bine diferentele dintre grafice am ales intervalul de timp care acopera numai doi ani, respectiv 2001-2002 in toate cele patru cazuri rezultatele obtinute referitor la semnul ratei de crestere a PiB-ului sunt apropiate, exceptie facand lunile din primul trimestru al fiecarui an in opinia noastra, cea mai indicata este metoda generala cu diferente constante luate zero in trimestrul i al fiecarui an, deoarece tine cont atat de evolutia seriilor lunare aditionale, cat si de influenta celorlati factori necunoscuti De asemenea, aceasta metoda minimizeaza efectul nefavorabil determinat de diferentele mari dintre trimestrul i si trimestrul iV al anului anterior, atribuit factorilor necunoscuti Astfel, seria PiB2 determinata cu aceasta metoda aproximeaza cel mai bine trendul PiB- ului lunar Exista cateva avantaje imediate ale metodei algebrice Comparativ cu metoda econometrica aplicata la seriile de timp in preturi comparabile, metoda algebrica estimeaza PiB-ul lunar direct in preturi curente si utilizeaza in algoritmul de interpolare un numar mai mare de serii aditionale De exemplu, seria lunara a cheltuielilor guvernamentale, de care ne-am folosit la interpolarea consumului guvernamental, nu mai este operabila pentru analiza econometrica; la fel se intampla cu taxele nete 27 Figura nr 3 6 Evolutia PiB-ului lunar estimat prin cele patru modele algebrice 200000 180000 160000 140000 120000 100000 80000 60000 40000 20000 0 123456789101112131415161718192021222324 Pe de alta parte, o metoda statistica nu da rezultate satisfacatoare la estimarea seriilor in valori absolute, asa cum am procedat in cazul metodei algebrice, datorita diferentelor mari intre valorile extreme ale seriei Aceste diferente acorda o pondere mai mare informatiei incluse in ultimele observatii atunci cand se estimeaza parametrii modelului Metoda Chow si Lin Seriile in preturi constante ale anului 1997 care se pot interpola prin metode econometrice sunt PiB-ul exclusiv agricultura, notat PiBA, si consumul final, notat CONSF Valoarea adaugata bruta din agricultura, notata VABA, se interpoleaza separat pe baza unor coeficienti de impartire a valorilor lunare pe trimestru, stabiliti in urma consultarii cu specialistii de la institutul de Economie Agrara Exportul net se determina, direct ca si in cazul metodei algebrice, pornind de la seriile lunare ale exportului si importului de bunuri si servicii din balanta de plati Pentru a transforma indicatorii de comert exterior (exprimati in mil s USA) in preturile medii ale anului 1997 vom aplica o rata de schimb calculata din seriile trimestriale (asa cum este de asteptat aceasta rata empirica fluctueaza in jurul unei valori fixe) PiB-ul lunar este suma elementelor PiBA si VABA, din care se obtine ulterior seria lunara a formarii brute de capital FBC prin scaderea consumului final si exportului net Conform statisticii sistemului conturilor nationale, seriile aditionale disponibile in preturi constante sunt cifra de afaceri a intreprinderilor din comertul cu amanuntul RS, potrivita pentru interpolarea consumului final, si respectiv PiB1 PiB2 PiB3 PiB4 Metoda directa cu variabile aditionale - PiB1 Metoda generala cu diferente constante luate zero in trimestrul i al fiecarui an - PiB2 Metoda generala cu diferente constante luate zero in trimestrul i al primului an - PiB3 Metoda directa fara variabile aditionale - PiB4 28 productia industriala PRODi si cifra de afaceri RS pentru interpolarea PiB-ului exclusiv agricultura Toate seriile se calculeaza in valori absolute prin inmultirea indicilor de volum cu valoarea nominala medie din anul 1997 Cheltuielile guvernamentale si taxele nete sunt eliminate din model deoarece nu avem un deflator pentru transformarea lor in preturi comparabile in figurile alaturate sunt reprezentate grafic evolutiile seriilor trimestriale PiBA si CONSF alaturi de PRODi si RS Pentru a se vedea mai bine corelatia seriei RS cu celelalte serii, am multiplicat RS cu 3 Figura nr 3 7 Evolutia PiB-ului exclusiv agricultura in corelatie cu productia industriala si cifra de afaceri a intreprinderilor din comertul cu amanuntul (multiplicata cu 3) 80000 70000 60000 50000 40000 30000 20000 10000 0 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 PiB-ul fara agricultura, Productia industriala si Cifra de afaceri in pretturi constante ale anului 1997 in perioada 1997 2003 Figura nr 3 8 Evolutia consumului final in corelatie cu cifra de afaceri a intreprinderilor din comertul cu amanuntul, multiplicata cu 3 90000 80000 70000 60000 50000 40000 30000 20000 10000 0 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 Q1 Q3 29 Din figurile nr 3 7 si 3 8 se observa existenta unei corelatii intre indicatori Analizele econometrice arata imbunatatirea rezultatelor in cazul ecuatiei de regresie a consumului final daca se considera productia industriala ca variabila explicativa alaturi de cifra de afaceri din comert Astfel, estimam ecuatiile de regresie trimestriale ale PiB-ului exclusiv agricultura (PiBA) si consumului final (CONSF) in functie de productia industriala si cifra de afaceri din comert, in ambele cazuri Totodata, din motive statistice, datele trimestriale ale anului 1997 induc erori in algoritmul de estimare al ecuatiilor de regresie, ceea ce ne determina sa restrangem perioada observatiilor care se folosesc la estimarea modelului la intervalul 1998:01 - 2003:04 Rezultatele regresiilor trimestriale sunt descrise in tabelul nr 2 si tabelul nr 3 Tabelul nr 2 Regresia PiB-ului exclusiv agricultura in functie de productia industriala si cifra de afaceri din comert Dependent Variable: PiBA Method: Least Squares Date: 06 15 04 Time: 06:26 Sample: 1998:1 - 2003:4 included observations: 24 Variable Coefficient Std Error t-Statistic Prob C -26122,17 7069,181 -3,695219 0,0013 PRODi 0,932916 0,120563 7,738002 0,0000 RS 2,595695 0,438367 5,921277 0,0000 R-squared 0,860805 Mean dependent var 53776,47 Adjusted R-squared 0,847548 S D dependent var 6244,416 S E of regression 2438,136 Akaike info criterion 18,55232 Sum squared resid 1,25E+08 Schwarz criterion 18,69958 Log likelihood -219,6279 F-statistic 64,93376 Durbin-Watson stat 0,963564 Prob(F-statistic) 0,000000 PiBA = -26122,17234 + 0,9329159531*PRODi + 2,595694741*RS Valorile mici ale parametrului Durbin Watson (0,96) indica prezenta corelatiei pozitive intre erori Valoarea coeficientului de determinare al regresiei de 0,86 arata ca seriile productiei industriale si cifrei de afaceri din comert explica satisfacator dinamica PiB-ului, exclusiv agricultura Acest rezultat este redat grafic in figura nr 3 9 unde se poate observa evolutia seriei reale si a seriei calculate de model, alaturi de graficul erorilor 30 Regresia trimestriala PiBA Figura nr 3 9 Tabelul nr 3 Regresia consumului final in functie de productia industriala si cifra de afaceri din comert Dependent Variable: CONSF Method: Least Squares Date: 06 15 04 Time: 07:23 Sample: 1998:1 - 2003:4 included observations: 24 Variable Coefficient Std Error t-Statistic Prob C -54462,64 14995,41 -3,631953 0,0016 PRODi 0,821787 0,255743 3,213338 0,0042 RS 5,410533 0,929882 5,818518 0,0000 R-squared 0,737096 Mean dependent var 57861,09 Adjusted R-squared 0,712058 S D dependent var 9638,165 S E of regression 5171,865 Akaike info criterion 20,05632 Sum squared resid 5,62E+08 Schwarz criterion 20,20358 Log likelihood -237,6759 F-statistic 29,43854 Durbin-Watson stat 1,487354 Prob(F-statistic) 0,000001 CONSF = -54462,64304 + 0,8217874332*PRODi + 5,410532634*RS 31 Aceasta regresie se imbunatateste din punct de vedere al comportamen- tului erorilor, asa cum o demonstreaza parametrul Durbin Watson (1,48) Valoarea coeficientului de determinare nu foarte ridicata (0,73) arata ca exista si alti factori importanti care explica evolutia consumului final CONSF Regresia trimestriala CONSF Figura nr 3 10 Estimarea seriilor lunare in ambele cazuri cu regresiile lunare avand coeficientii si erorile calculate din regresiile trimestriale este imbunatatita de conditia insumarii valorilor lunare la valoarea trimestriala (metoda Chow si Lin ia in calcul aceasta conditie) Astfel, in cazul seriei lunare PiBA estimate cu regresia Chow si Lin cu parametrii determinati din regresia cu variabile trimestriale, am obtinut un procent al erorilor fata de valorile estimate situat in intervalul [0 3%, 6,5%] pentru perioada 1998-2003 in ceea ce priveste seria lunara CONSF procentul erorilor fata de valorile estimate sunt mai mari, atingand valori procentuale de 10%-12% in trimestrele i si iV interpolarea PiB-ului exclusiv agricultura este mai precisa decat cea a consumului final; de asemenea precizia maxima a rezultatelor se obtine daca aplicam metoda econometrica pentru intervalul trimi 1999-trimiV 2003 Precizam ca rezultatele modelului mai putin multumitoare pentru prima perioada a anului 1997 se datoreaza schimbarii metodologiei statistice de evaluare a seriilor lunare PiB-ul lunar este suma dintre VABA (stabilita in afara modelului) si PiBA (a se vedea graficul din figura nr 3 11) Seria este stationara si are un caracter sezonier pronuntat PiB-ul a fost construit pe baza insumarii influentelor lunare ale productiei agricole, productiei industriale si cifrei de afaceri a intreprinderilor din comertul cu amanuntul Evolutia la nivel lunar este modulata de tendintele medii cunoscute ale seriilor trimestriale 32 Figura nr 3 11 Evolutia PiB-ului lunar in preturile medii ale anului 1997 in perioada 1997-2003 35000 30000 25000 20000 15000 10000 5000 0 PiB-ul lunar in preturi constante in perioada 1997-2003 Valorile maxime ale PiB-ului lunar in preturi constante sunt atinse in lunile iulie si octombrie (in cazul PiB-ului lunar in preturi curente valoarea maxima este atinsa in luna octombrie, dar un salt rapid se realizeaza si in luna iulie) Valorile minime sunt in lunile martie-1997, februarie-1998, ianuarie-ceilalti ani in figura nr 3 12 am reprezentat grafic evolutia PiB-ului lunar fara agricultura, PiBA, in perioada 1997-2003 Evolutia acestui indicator este mai putin oscilanta decat cea a PiB-ului Din comparatia celor doi indicatori VABA (valoarea adaugata bruta din agricultura) si PiBA constatam ca valorile minime ale PiB-ului sunt determinate de PiBA in timp ce valorile maxime sunt determinate de VABA pentru luna iulie si de VABA, PiBA pentru luna octombrie Scenariile diferite pe care le-am aplicat privind evolutia lunara a valorii adaugate brute din agricultura nu afecteaza rezultatul de maxim al PiB-ului in luna octombrie Figura nr 3 12 Evolutia PiB-ului lunar (fara agricultura) in preturile medii ale anului 1997, in perioada 1997-2003 25000 20000 15000 10000 5000 0 33 Modelul de prognoza al PiB-ului lunar in economia romaneasca Seria lunara a PiB-ului estimata cu metoda Chow si Lin pentru perioada 1997 - 2003 (Figura nr 3 11) se poate utiliza pentru construirea unui model care sa prognozeze PiB-ul lunar pe o perioada de un trimestru (trimestrul i, 2004) in care sunt disponibile datele lunare ale seriilor aditionale (productia industriala si cifra de afaceri din comertul cu amanuntul) in acest mod avem posibilitatea sa estimam valorile lunare ale PiB-ului si valoarea trimestriala cu trei-patru luni inainte de publicarea datelor oficiale de catre institutul National de Statistica Deoarece datele trimestriale ale anului 2003 sunt provizorii, nu este indicat sa facem prognoze pe o perioada mai lunga de un trimestru datorita propagarii erorilor continute in datele provizorii oficiale De asemenea, pentru a cerceta stabilitatea modelului, reducem seria istorica a PiB-ului cu trei luni, reestimam modelul folosind datele lunare partiale pana in trimestrul iii 2003 (inclusiv) si prognozam PiB-ul pe lunile trimestrului iV din anul 2003 Valorile PiB-ului lunar cunoscute din trimestrul iV 2003 (determinate cu metoda Chow si Lin) si valorile prognozate din trimestrul i 2004 folosind modelul estimat cu date lunare complete trebuie sa se pozitioneze in intervalul de plus minus o abatere standard calculat de modelul estimat cu datele lunare partiale Pentru a construi modelul de prognoza al PiB-ului lunar optam, ca si in cazul seriilor PiB-ului exclusiv agricultura si a consumului final, pe corelatia dintre dinamica PiB-ului si dinamica seriilor lunare ale productiei industriale si cifrei de afaceri din comertul cu amanuntul Astfel, vom estima o ecuatie de regresie cu date lunare avand forma generala PiBt=c1+c * PRODit+ c * RSt + F + ε t t12t3t tt unde momentul t ia valori intre 1 si T (nu cunoastem inca dimensiunea probei statistice pe care o folosim la estimarea regresiei), ε t este termenul eroare, c1,c2, c3 sunt coeficientii de estimat Functia F t este introdusa pentru a lua in calcul efectele induse de sezonalitatea lunara Prima dificultate legata de specificarea modelului vizeaza alegerea intervalului statistic optim care sa conduca la regresia cea mai buna pentru realizarea prognozelor Deoarece nu avem serii observabile decat la nivel trimestrial, am realizat mai multe tipuri de modele ale PiB-ului trimestrial rezultand ca cea mai buna regresie este aceea estimata in intervalul [trim i, 1999 - trim iV, 2003] Se constata faptul ca anii 1997 si 1998 inrautatesc corelatia dintre seriile PiB-ului, productiei industriale si cifrei de afaceri din comert, si din acest motiv ii excludem din proba statistica atunci cand estimam regresia lunara Rezultatele regresiei trimestriale sunt descrise in tabelul nr 4 si figura nr 3 13 Mentionam ca desi exista o sezonalitate pronuntata a PiB-ului trimestrial datorita agriculturii, rezultatele regresiei sunt pozitive, o parte din sezonalitatea PiB-ului fiind preluata de sezonalitatea celor doua serii PRODi si RS 34 Tabelul nr 4 Regresia PiB-ului trimestrial in functie de productia industriala si cifra de afaceri din comert Dependent Variable: PiB Method: Least Squares Date: 06 16 04 Time: 17:16 Sample: 1999:1 2003:4 included observations: 20 Variable Coefficient Std Error t-Statistic Prob C -70615,03 13876,99 -5,088642 0,0001 PRODi 0,794277 0,260797 3,045572 0,0073 RS 7,269429 0,963141 7,547626 0,0000 R-squared 0,858058 Mean dependent var 64748,84 Adjusted R-squared 0,841359 S,D, dependent var 11879,07 S E of regression 4731,404 Akaike info criterion 19,89931 Sum squared resid 3,81E+08 Schwarz criterion 20,04867 Log likelihood -195,9931 F-statistic 51,38370 Durbin-Watson stat 1,256958 Prob(F-statistic) 0,000000 Figura nr 3 13 Regresia PiB-ului trimestrial in perioada [trim i, 1999 - trim iV, 2003] 35 Cea de-a doua dificultate este data de sezonalitatea seriilor la nivel lunar, evidentiata de noi in model prin introducerea functiei F Analiza mai multor tipuri t de modele ne-a condus la alegerea F= PiB Rezultatele regresiei lunare a tt-12 PiB-ului sunt descrise in tabelul nr 4 si figura nr 3 14 Acesta este si modelul pe care il vom folosi la realizarea de prognoze intervalul de estimare este t = [lun i 1999 - lun Xii 2003] 2 Valorile AiC, BiC si adjusted R (se observa o valoare foarte ridicata a lui 2 R ) sunt cele mai bune in comparat ie cu valorile obt inute pentru celelalte modele analizate Corelograma reziduurilor indica lipsa corelatiei Testele Durbin-Watson (pentru autocorelarea reziduurilor) si Goldfend-Quandt (pentru heteroskedasti- citate) indica la randul lor lipsa autocorelarii si a heteroskedasticitatii reziduurilor Din analizarea testelor statistice concluzionam ca valorile estimate pentru parametrii regresiei sunt semnificative si modelul este specificat corect Tabelul nr 5 Regresia PiB-ului lunar in functie de productia industriala si cifra de afaceri din comert Dependent Variable: PiB Method: Least Squares Date: 06 17 04 Time: 06:49 Sample: 1999:01 - 2003:12 included observations: 60 Variable Coefficient Std Error t-Statistic Prob C PRODi RS PiB(-12) -6011,468 0,341226 0,651184 0,917360 1206,150 -4,984012 0,061812 5,520432 0,293781 2,216560 0,037138 24,70157 0,0000 0,0000 0,0307 0,0000 R-squared 0,978251 Mean dependent var 21582,95 Adjusted R-squared 0,977086 S D dependent var 4302,361 S E of regression 651,2587 Akaike info criterion 15,86003 Sum squared resid 23751725 Schwarz criterion 15,99965 Log likelihood -471,8009 F-statistic 839,6300 Durbin-Watson stat 1,901972 Prob(F-statistic) 0,000000 PiB = -6011,46823 + 0,917360016*PiB(-12) + 0,3412261873*PRODi + 0,6511841537*RS 36 Regresia PiB-ului lunar in perioada [lun i, 1999 - lun Xii, 2003] Figura nr 3 14 Pe baza ecuatiei de regresie descrisa in tabelul nr 5 putem prognoza PiB-ul lunar in trimestrul i, 2004 Au fost gasite urmatoarele valori: Luna Simbol Valoare ianuarie 2004 PiB(61) 17937 97 Februarie 2004 PiB(62) 18104 24 Martie 2004 PiB(63) 21764 12 Produsul intern brut prognozat pentru perioada 1 i -31 iii 2004 este 57806,3 miliarde lei in preturi medii ale anului 1997, in crestere cu 8,7% fata de perioada 1 i -31 iii 2003 Datele oficiale ale institutului National de Statistica indica o crestere de 6,1% fata de acelasi trimestru din anul anterior Pentru a cerceta stabilitatea modelului la modificarea dimensiunii probei statistice, reducem seria istorica a PiB-ului cu trei luni si reestimam regresia folosind datele lunare din perioada [lun i 1999-lun iX 2003] Cu regresia obtinuta in acest mod prognozam PiB-ul pe lunile trimestrului iV din anul 2003, notandu-le cu PiB iV 03 Rezultatele sunt prezentate in figura nr 3 15 Pe acelasi grafic sunt prezentate valorile PiB-ului lunar din trimestrul i, 2004 prognozate cu regresia estimata cu proba statistica completa, notate cu PiB i 04 37 Figura nr 3 15 Regresia PiB-ului lunar: prognoza trim iii 2003 (proba partiala)-PiB iV 03 prognoza trim i 2004 (proba completa)-PiB i 04; si intervalul de incredere de plus minus o abatere standard al regresiei cu proba partiala 36000 32000 28000 24000 20000 16000 2003:07 2003:10 2004:01 Valorile cunoscute (efective) ale PiB-ului lunar obtinute cu metoda Chow si Lin pentru trimestrul iV 2003 (seria PiB) impreuna cu valorile lunare din trimes- trul i 2004 prognozate cu regresia estimata cu proba statistica completa (seria PiB i 04) se pozit ioneaza in intervalul de plus minus cu o abatere standard calculata cu regresia estimata cu proba partiala [lun i, 1999-lun iX, 2003] Aceasta situatie demonstreaza precizia prognozelor si calitatea ecuatiei de regresie (parametrii ecuatiei de regresie isi pastreaza proprietatile chiar daca au fost estimati utilizand probe statistice diferite) Concluzii • Modelul PiB-ului lunar este important pentru ca prezinta dinamica pe termen scurt a outputului ‘global' al economiei Astfel, se pot realiza principalele corelatii dintre indicatorii monetari-financiari si indicatorii macroeconomici, in scopul implementarii corecte a politicilor monetar-fiscale pe termen scurt tintirea inflatiei pe termen scurt, politica cursului de schimb si a ratei dobanzii nu pot fi realizate fara detinerea de informatii corecte asupra cresterii economice Rezultatul cresterii economice la nivel anual ar trebui vazut ca o rezultanta a comportamentelor decidentilor politici, pe de o parte, si a raspunsurilor economiei, pe de alta parte, ca o reactie si o adaptare continua la schimbarile de politici PiB SUP PiBiV03 iNF PiB i 04 38 • Modelul PiB-ului lunar genereaza o baza de date cu informat ii condensate Deoarece seria PiB-ului lunar prezinta un numar mare de observatii ea poate fi analizata cu ajutorul tehnicilor moderne econometrice ARiMA, ARCH, GARCH pentru a se desprinde astfel regularitatile in dinamica pe termen scurt de-a lungul unui an Aceste regularitati sunt specifice unui anumit tip de economii si nu se disting clar decat la frecvente inalte ale indicatorilor economici Compararea diverselor tipuri de economii de tranzitie scoate in evidenta deosebiri si asemanari utile pe baza indicatorilor lunari • Prin seria lunara a PiB-ului punem la dispozit ia utilizatorului un set complet de informatii privind evolutia PiB-ului potential De exemplu, PiB-ul potential ar trebui calculat (cu filtrul Hodrick Prescott) si la nivel trimestrial, iar acest lucru este posibil numai daca folosim seria lunara • Seria PiB-ului lunar se utilizeaza pentru a prognoza PiB-ul trimestrial pe un singur trimestru sau cel mult pe doua trimestre Aceste prognoze se bazeaza pe datele statistice publicate ale seriilor lunare precum productia industriala si cifra de afaceri a intreprinderilor din comert Deoarece PiB-ul trimestrial se publica cu o intarziere de trei-patru luni comparativ cu ceilalti indicatori lunari, avem astfel posibilitatea sa cunoastem valoarea estimativa a acestuia cat de curand posibil in studiul de fata noi am realizat o prognoza a PiB-ului pentru trimestrul i 2004 avand la dispozitie datele statistice lunare pana in luna aprilie • Scopul principal al studiului este estimarea seriei lunare a PiB-ului care interpoleaza datele trimestriale istorice in perioada 1997-2004, pentru a putea fi utilizata de beneficiar in modelele econometrice pe termen scurt Rezultatele obtinute cu metoda algebrica si metoda Chow si Lin arata clar ca precizia estimarilor depinde in mare masura de folosirea seriilor aditionale lunare corelate ipotetic (conform teoriei statisticii economice) cu seria care se interpoleaza Astfel, estimarile pe baza metodei algebrice au un grad de precizie mai ridicat, deoarece aceasta metoda are la indemana (spre deosebire de metoda econometrica Chow&Lin) cheltuielile guvernamentale si taxele nete lunare Conform estimarilor noastre pentru economia Romaniei, valoarea maxima a PiB-ului lunar, atat in preturi curente cat si in preturi constante, este atinsa in luna octombrie a fiecarui an, iar valoarea minima in luna ianuarie, cu unele exceptii: februarie (preturi constante-anul 1998) si martie (preturi constante-anul 1997) Seriile PiB-ului lunar reprezentative estimate prin metoda algebrica (metoda generala cu diferente constante luate zero in trimestrul i al fiecarui an) si metoda econometrica sunt incluse in anexele nr 1 si 2 39 Anexa nr 1 Seria reprezentativa a PiB-ului lunar in preturi curente estimata prin metoda algebrica (metoda generala cu diferente constante luate zero in trimestrul i al fiecarui an): perioada PiB industrie Agricultura Celelalte impozite 1997 ian 11797 3726,3 553,2 6084,9 1432,3 Feb 12590 4319,2 632,2 6215,9 1422,6 Mar 15394 6641,6 790,3 6400,6 1562,0 Apr 15988 6295,6 831,6 7374,5 1486,3 Mai 16974 6523,9 942,5 7907,8 1599,5 iun 20407 6427,6 3769,9 8603,5 1605,6 iul 25117 6583,5 7594,4 9255,2 1683,6 Aug 24125 6624,0 6036,6 9926,1 1538,4 Sep 25224 6634,2 5841,9 10638,8 2108,9 Oct 30721 8022,3 9270,2 11458,2 1970,3 Nov 28241 8713,3 5562,1 12112,2 1853,8 Dec 26348 7582,1 3708,1 13213,4 1844,7 1998 ian 20697 6914,4 720,2 10783,4 2279,1 Feb 21076 7157,9 823,0 10761,7 2333,3 Mar 22904 8028,6 1028,8 10983,1 2863,3 Apr 23838 7893,9 1034,6 12438,6 2470,7 Mai 25934 8695,6 1172,5 13411,0 2655,0 iun 31620 9309,4 4690,0 14652,0 2968,9 iul 35779 8788,6 8652,9 15486,0 2851,1 Aug 34642 8508,4 6877,9 16302,0 2953,9 Sep 37037 9278,2 6656,0 17139,7 3963,3 Oct 42721 9833,9 11054,0 17765,4 4067,3 Nov 38073 9521,6 6632,4 18415,1 3504,0 Dec 36873 9123,1 4421,6 19650,6 3677,8 1999 ian 28065 7343,6 972,5 16275,7 3473,2 Feb 28697 7808,5 1111,4 16289,3 3487,5 Mar 32516 9356,9 1389,3 16614,1 5155,5 Apr 34088 9637,6 1357,0 18382,9 4710,2 Mai 36680 10589,4 1538,0 19711,6 4841,2 iun 44092 11613,3 6151,9 21182,0 5145,2 iul 51486 11990,4 10988,6 22799,8 5707,1 Aug 51806 12442,8 8734,5 24282,7 6345,9 Sep 54761 13515,2 8452,8 25750,9 7042,4 40 perioada PiB industrie Agricultura Celelalte impozite Oct 63774 13499,9 15521,0 28075,0 6677,8 Nov 60790 14335,5 9312,6 30511,4 6630,7 Dec 58975 13210,6 6208,4 33517,8 6038,4 2000 ian 41769 11296,1 1253,4 23578,2 5641,3 Feb 43298 12845,3 1432,5 23454,7 5565,1 Mar 47229 14977,1 1790,6 24118,0 6343,9 Apr 50801 14967,8 2073,4 27298,6 6461,1 Mai 55374 17321,2 2349,9 28980,7 6722,4 iun 67562 19392,7 9399,5 31394,4 7375,1 iul 74840 19680,9 13953,1 34036,1 7169,5 Aug 75983 19289,4 11090,9 36806,4 8796,5 Sep 80284 19592,4 10733,1 40036,6 9921,8 Oct 90887 21793,6 16924,3 42032,6 10136,3 Nov 89489 24866,8 10154,6 44451,6 10016,1 Dec 86258 23456,4 6769,7 47732,7 8298,8 2001 ian 63621 19141,0 1802,8 35703,5 6973,7 Feb 64004 20660,7 2060,4 35011,2 6272,0 Mar 69248 23512,9 2575,5 35610,7 7549,3 Apr 76515 24419,5 3300,8 40481,5 8313,5 Mai 84280 28397,2 3740,9 43208,9 8933,0 iun 100687 28810,4 14963,8 47031,6 9881,0 iul 111109 28986,7 24607,3 46722,5 10792,2 Aug 105948 28287,6 19559,7 46486,2 11614,6 Sep 105468 27928,9 18928,7 46448,3 12162,4 Oct 130618 32160,3 30836,1 53552,5 14069,1 Nov 127973 35106,8 18501,7 60960,4 13404,5 Dec 127770 31515,0 12334,5 69306,2 14614,8 2002 ian 85248 26853,4 2232,8 46676,1 9485,9 Feb 85963 28443,3 2551,8 45344,9 9623,1 Mar 89565 30626,4 3189,8 46262,0 9486,8 Apr 100067 33367,0 4090,8 51122,6 11486,6 Mai 107746 36041,7 4636,2 54264,4 12804,0 iun 128012 38247,9 18544,7 57894,0 13325,7 iul 139902 39553,2 25647,6 59708,4 14992,8 Aug 134352 37737,7 20386,5 61417,0 14810,3 Sep 136723 38583,1 19728,9 63083,2 15327,6 Oct 169076 43593,7 34852,6 74080,6 16548,6 41 perioada PiB industrie Agricultura Celelalte impozite Nov 166876 45741,8 20911,6 85111,1 15111,4 Dec 168727 40701,5 13941,0 97668,4 16416,1 2003 ian 106980 34089,6 2525,0 58027,8 12338,0 Feb 106651 35113,2 2885,7 56323,6 12328,3 Mar 114072 39500,0 3607,1 57080,2 13885,1 Apr 124379 41113,0 4813,6 62941,7 15511,2 Mai 132009 44114,8 5455,4 66459,8 15979,4 iun 153268 46202,8 21821,6 70695,3 14548,0 iul 180326 49393,6 38344,7 73362,0 19226,2 Aug 166994 45956,0 30479,1 75636,3 14922,3 Sep 179768 51828,1 29495,9 78096,3 20347,5 Oct 210225 53674,2 40841,2 94615,8 21094,2 Nov 206988 51236,2 24504,7 111258,6 19988,9 Dec 209116 44952,0 16336,5 130253,1 17574,8 PiB Consum Gospodarii Adm Publica Adm Privata FBC EXP net 1997 ian 11797 10652,4 9361,5 1214,7 76,2 2064,1 -919,9 Feb 12590 11737,5 10210,9 1450,5 76,2 1176,1 -323,6 Mar 15394 13901,1 11330,7 2494,2 76,2 1779,0 -285,7 Apr 15988 14902,6 12637,5 2177,6 87,5 2057,1 -971,6 Mai 16974 16277,3 13946,9 2243,0 87,5 2111,1 -1414,7 iun 20407 17743,7 15126,3 2529,9 87,5 4245,6 -1582,7 iul 25117 19150,4 16907,8 2139,6 103,1 7974,8 -2008,5 Aug 24125 19698,5 17293,3 2302,1 103,1 4450,2 -23,6 Sep 25224 22861,4 19138,2 3620,2 103,1 3246,3 -884,0 Oct 30721 23412,6 20070,5 3148,2 193,8 9162,4 -1854,1 Nov 28241 22217,1 19278,6 2744,6 193,8 7425,8 -1401,5 Dec 26348 26065,0 20936,1 4935,1 193,8 6478,8 -6195,6 1998 ian 20697 18977,7 16377,9 2397,2 202,7 2377,4 -658,0 Feb 21076 18921,1 15627,5 3090,9 202,7 2410,4 -255,5 Mar 22904 21462,3 17842,7 3416,9 202,7 3051,2 -1609,6 Apr 23838 22679,4 19306,0 3163,7 209,7 3300,3 -2141,9 Mai 25934 24861,1 20749,1 3902,3 209,7 3674,6 -2601,5 iun 31620 27528,6 22556,3 4762,5 209,7 6279,9 -2188,3 iul 35779 26779,2 24413,5 2099,0 266,7 11499,5 -2500,1 Aug 34642 30400,9 25824,7 4309,5 266,7 5510,9 -1269,7 Sep 37037 36461,8 28975,6 7219,5 266,7 4174,0 -3598,5 Oct 42721 36326,9 29477,7 6426,5 422,7 9794,6 -3400,9 Nov 38073 33478,7 27870,7 5185,3 422,7 7824,4 -3230,0 42 PiB Consum Gospodarii Adm Publica Adm Privata FBC EXP net Dec 36873 36794,7 29604,5 6767,6 422,7 6628,0 -6549,7 1999 ian 28065 25740,9 21476,6 3961,7 302,5 3638,2 -1313,9 Feb 28697 26978,1 21851,0 4824,5 302,5 2342,4 -623,7 Mar 32516 33200,2 27372,5 5525,2 302,5 1028,6 -1712,9 Apr 34088 33211,1 27675,8 5244,2 291,1 3751,4 -2874,7 Mai 36680 34968,8 28997,1 5680,6 291,1 3929,4 -2218,1 iun 44092 38073,5 30172,0 7610,4 291,1 8369,5 -2350,6 iul 51486 42298,2 34577,6 7239,3 481,3 10405,2 -1217,5 Aug 51806 42033,5 35895,3 5656,9 481,3 10430,4 -658,0 Sep 54761 46877,6 39595,5 6800,7 481,3 8851,0 -967,3 Oct 63774 49532,1 40922,1 7672,8 937,2 15822,6 -1580,9 Nov 60790 52049,8 42811,1 8301,5 937,2 12854,3 -4113,9 Dec 58975 59397,7 47938,5 10521,9 937,2 6317,7 -6740,2 2000 ian 41769 34965,0 27685,2 6902,3 377,5 7421,7 -617,7 Feb 43298 37747,8 30780,8 6589,4 377,5 5131,4 418,5 Mar 47229 43002,2 35027,7 7597,0 377,5 6522,2 -2295,0 Apr 50801 47284,3 38348,2 8511,2 424,8 5700,7 -2184,0 Mai 55374 49371,8 40375,0 8572,0 424,8 13125,9 -7123,6 iun 67562 54180,4 42643,6 11112,0 424,8 15141,7 -1760,5 iul 74840 57886,8 48071,5 9183,2 632,1 19757,0 -2804,2 Aug 75983 62014,3 50740,0 10642,2 632,1 15274,6 -1305,7 Sep 80284 68706,4 55491,1 12583,2 632,1 13385,9 -1808,4 Oct 90887 74904,0 57639,9 15596,1 1668,1 22130,8 -6148,1 Nov 89489 75228,6 59920,8 13639,7 1668,1 22919,8 -8659,3 Dec 86258 87241,1 67150,4 18422,7 1668,1 9979,2 -10962,7 2001 ian 63621 55757,9 44295,1 10905,5 557,2 12389,1 -4525,9 Feb 64004 56854,5 46046,3 10250,9 557,2 14534,5 -7384,7 Mar 69248 64945,6 53628,7 10759,6 557,2 8141,3 -3838,6 Apr 76515 69617,6 58517,3 10518,2 582,2 18086,3 -11188,6 Mai 84280 76488,2 63342,8 12563,2 582,2 18965,2 -11173,4 iun 100687 83947,6 64056,2 19309,2 582,2 21369,0 -4629,9 iul 111109 96315,7 68746,3 26791,9 777,5 20207,3 -5414,3 Aug 105948 76389,4 65180,6 10431,4 777,5 32403,3 -2844,6 Sep 105468 75404,6 64251,0 10376,1 777,5 29641,8 421,9 Oct 130618 94809,8 77624,2 15014,9 2170,7 49055,9 -13247,7 Nov 127973 108699,5 88395,7 18133,2 2170,7 29120,0 -9846,2 Dec 127770 134976,0 104048,7 28756,7 2170,7 9614,2 -16819,8 2002 ian 85248 73312,7 58885,0 13721,0 706,7 17577,8 -5642,3 Feb 85963 76065,5 62196,5 13162,4 706,7 15765,0 -5867,5 Mar 89565 85135,2 70884,8 13543,7 706,7 9968,5 -5538,8 Apr 100067 97232,5 81153,5 15401,7 677,3 12135,5 -9300,9 Mai 107746 91606,9 76045,4 14884,2 677,3 23745,1 -7605,8 43 PiB Consum Gospodarii Adm Publica Adm Privata FBC EXP net iun iul 128012 139902 96480,2 107204,6 75117,5 85648,6 20685,4 20570,9 677,3 985,2 37985,2 41895,5 -6453,1 -9198,2 Aug 134352 99046,2 81599,9 16461,1 985,2 36313,5 -1008,2 Sep 136723 106163,0 81856,2 23321,6 985,2 33925,3 -3365,6 Oct 169076 120728,1 95056,7 23168,7 2502,7 59263,4 -10915,9 Nov 166876 134215,7 105263,0 26450,0 2502,7 43026,1 -10365,8 Dec 168727 163142,9 126028,9 34611,3 2502,7 17862,0 -12277,9 2003 ian 106980 90448,4 77186,4 12493,9 768,0 20167,0 -3635,1 Feb 106651 93670,3 77111,0 15791,3 768,0 16927,0 -3946,5 Mar 114072 103160,9 85969,7 16423,3 768,0 17920,2 -7008,8 Apr 124379 115335,4 97199,9 17371,8 763,7 23230,6 -14186,5 Mai 132009 116222,9 98939,9 16519,3 763,7 30652,4 -14866,0 iun 153268 123135,0 100039,4 22331,8 763,7 42786,5 -12653,8 iul 180326 132185,0 108553,0 22622,6 1009,3 57720,8 -9579,3 Aug 166994 131765,5 111050,9 19705,3 1009,3 40997,0 -5768,7 Sep 179768 142476,9 113255,9 28211,7 1009,3 51342,7 -14051,8 Oct 210225 160277,7 129336,9 28462,3 2478,5 73449,3 -23501,7 Nov 206988 165866,9 132004,6 31383,7 2478,5 59610,4 -18488,9 Dec 209116 199991,3 157607,2 39905,6 2478,5 31175,5 -22050,4 Sursa: institutul National de Statistica si calcule proprii Asistenta din partea specialistilor institutului National de Statistica: Silvia Caragea; Gheorghe Mihai 44 Anexa 2 Seria reprezentativa a PiB-ului lunar in preturi constante ale anului 1997 estimata prin metoda econometrica Chow&Lin perioada PiB PiB exclusiv Agric Agricultura 1997 ian 18492,5 17704,4 788,1 Feb 17663,3 16762,5 900,7 Mar 16334,8 15208,9 1125,9 Apr 18027,0 17044,1 982,9 Mai 18245,9 17132,0 1113,9 iun 21192,3 16736,6 4455,7 iul 24675,7 17301,8 7374,0 Aug 22919,9 17058,5 5861,4 Sep 23236,6 17564,3 5672,3 Oct 27801,9 19172,9 8629,0 Nov 24148,9 18971,5 5177,4 Dec 20186,9 16735,3 3451,6 1998 ian 15580,4 14906,0 674,4 Feb 14447,8 13677,0 770,7 Mar 17147,7 16184,3 963,4 Apr 17141,0 16239,5 901,6 Mai 17465,9 16444,1 1021,8 iun 20831,2 16744,1 4087,0 iul 22987,0 16379,5 6607,6 Aug 21617,7 16365,5 5252,2 Sep 23548,4 18465,7 5082,7 Oct 27016,0 19292,7 7723,3 Nov 22339,8 17705,8 4634,0 Dec 20618,3 17529,0 3089,3 1999 ian 14128,7 13510,5 618,2 Feb 14393,4 13686,9 706,5 Mar 18049,5 17166,3 883,2 Apr 17122,9 16229,3 893,6 Mai 17396,9 16384,1 1012,8 iun 20333,5 16282,4 4051,2 iul 23728,9 16963,6 6765,3 Aug 21948,4 16570,8 5377,5 Sep 22485,4 17281,4 5204,1 Oct 25815,2 17496,2 8319,1 Nov 23179,2 18187,8 4991,4 Dec 21042,7 17715,0 3327,6 2000 ian 13327,3 12752,5 574,8 Feb 15411,4 14754,5 656,9 Mar 18487,7 17666,5 821,1 Apr 16749,8 15832,5 917,3 Mai 18072,6 17032,9 1039,6 45 perioada PiB PiB exclusiv Agric Agricultura iun 21676,7 17518,2 4158,5 iul 23839,9 18435,0 5404,9 Aug 22497,9 18201,7 4296,2 Sep 22826,4 18668,8 4157,6 Oct 25942,1 19692,2 6249,9 Nov 23957,2 20207,3 3749,9 Dec 21606,6 19106,7 2500,0 2001 ian 14874,3 14293,7 580,6 Feb 15682,3 15018,7 663,6 Mar 18793,8 17964,3 829,5 Apr 18083,6 17107,2 976,4 Mai 19642,6 18536,0 1106,6 iun 21704,9 17278,6 4426,3 iul 25820,3 18886,7 6933,6 Aug 23833,2 18321,9 5511,3 Sep 23817,4 18483,8 5333,5 Oct 28672,2 20237,9 8434,4 Nov 25395,2 20334,6 5060,6 Dec 21470,8 18097,0 3373,7 2002 ian 15484,1 14881,6 602,4 Feb 16558,8 15870,4 688,5 Mar 18875,3 18014,7 860,6 Apr 20008,5 18936,5 1072,0 Mai 19564,6 18349,6 1214,9 iun 23184,9 18325,3 4859,7 iul 27026,9 20752,1 6274,8 Aug 24603,0 19615,4 4987,7 Sep 25044,2 20217,5 4826,8 Oct 30221,3 22134,8 8086,5 Nov 26439,0 21587,2 4851,9 Dec 22957,9 19723,3 3234,6 2003 ian 16754,2 16109,0 645,2 Feb 16814,1 16076,7 737,4 Mar 19590,3 18668,6 921,7 Apr 20113,4 19071,4 1042,0 Mai 20785,1 19604,2 1180,9 iun 24495,3 19771,7 4723,6 iul 28588,5 21880,7 6707,8 Aug 25674,3 20342,5 5331,8 Sep 26628,5 21468,7 5159,8 Oct 31773,6 23595,1 8178,5 Nov 26994,6 22087,5 4907,1 Dec 24985,2 21713,8 3271,4 Sursa: institutul National de Statistica si calcule proprii 46 Bibliografie Bernanke B ,Gertler M , Watson M - "Systemic monetary policy and the effects of oil price shocks", Brookings Papers on Economic Activity,1997 Chow G and Lin l - "Best linear unbiased interpolation, distribution, and extrapolation of time series by related series", The Review of Economics and Statistics, 1971 Chow G and Lin l - "Best linear unbiased estimation of missing observation in an economic time series", Journal of the American Statistical Association, September 1976 Harvey, Andrew C - Forecasting, Structural Time Series Models and the Kalman filter; Cambridge University Press, 1994 Hamilton, James D - Time Series Analysis, Princeton University Press, 1994 Kalman R - "A new approach to linear Filtering and prediction problems", Journal of Basic Engineering, Transactions of the ASME Series D,1960 Lanning S - "Missing observations: A simultaneous approach versus interpolation by related series", Journal of Economic and Social Measurement, April 1986 Leamer E - Model choice and specification analysis Handbook of Econometrics i, ed by Z Griliches and M D intriligator, North-Holland, 1983 Mills, Terence C - The Econometric Modelling of Financial Time Series; Cambridge University Press, 1999 